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FACILITY: MATH LIBRARY 


"ABSTRACT: 
This module contains routine MTH$CGSQRT - compute 
G COMPLEX*16 square root. 


VERSION: 1 
HISTORY: 


AUTHOR: 
Steven B. Lionel, 24-July-1979 


MODIFIED BY: 


SOOOGOSOSOOOSOOSOOOSOSOSOSOOSOSSOOOSOOSOSOSOSOSSOOOOSOOOOOOOOOOOOOO 
SOOOSOOSOOOSOSOOSOSOSOSSSSSSOSOSOOSOOOSOSOSOSOSOSOOOOOOOOSOOSOOSoOOo 
SOSOOSSOSSOSSOSSOSOSOSOSSOSSSSOSSOSOSOSOSOSOSOSOSOOOSOSOSOOOOOOOOoOOo 
SOOCOCOOSOSOOGSSOOOSSOOSSSSSSOOOCOSOOOCCOCOOOOOOCOOCOOOOOOoO 
NAME WN SO OD NAME WIN OOO NAUE WWI 2 O OONO UE WS OVODONOU ES Win 
Sete Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge 


BPS SF FWWIWWIWIMA WINING PNUPINPINUNPUNINIP 2 OO OO OO 


mM 5 , 
= W1sTORY 5 detailed Current att nistory TECSEDSISRE T:SRETY YAMA. SREGRIHESSBBr mansn "29% Gy) | 
6 -SBTTL HISTORY ; Detailed Current Edit History 


; Edit History 


: 1-001 - Adapted from MTHSCSQRT version ing f ve gy erties 
; 1-002 = Use general mode addressing. SBL 30-Nov-1981 
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»~SBTTL DECLARATIONS 
; INCLUDE FILES: 

: EXTERNAL pares: 

LEXTRN MTHSCGABS 
-EXTRN = MTHSGSQRT_RS 


; MACROS: 
NONE 


OOOOOOCOCOOOOCOOCOooO 


PSECT eet ee ATIONS: 
-PSECT _MTHSCODE PIC, SHR, LONG, EXE, NOWRT 


0000 


; EQUATED SYMBOLS: 


CSOOOCOSCOCSOOOCOCOOCOOOOOOOOOOOOOO 


COO O9 NINN NI NINN FAO OOO AOOOOUIUTIN 


OSOoooooooooooooo 
MO OONAOUE WN" OOONOAUS WOON 
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SOoOoooooooooooooooo 


SoOoooooooooo 


> OWN STORAGE: 
NONE 


hee £10: MT 
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0 ¢ ~SBTTL MTHSCGSQRT = compute G COMPLEX*16 square root 
399 
$ : FUNCTIONAL DESCRIPTION: 
§ : The square root of a complex number (r, i) is computed 
oy ; as follows: 
0 35 : ROOT = SQRT((ABS(r) + CABS((r, 1))) 7 2) 
3) 93 ; Q= i / (2*ROOT) 
94; 
$9 3 ; r i CSQRT((r, i)) 
000 98 ; 
000 99 ; >=0 an (ROOT, Q) 
000 100: <0 56 (Q, ROOT) 
000 101 : <0 <0 (-@, -ROOT) 
B88 1 ¢ 3 
444 HF ; CALLING SEQUENCE: 
9000 195 : CALL MTHSCGSQRT (result.wgc.r, arg.rgc.r) 
0000 109 : INPUT PARAMETERS: 
00000008 9009 109° 8 The G COMPLEX#16 t d 
ar = : * e 
0000 110 . 3 by veterence. asihaahn ote 
0000 111; 
0000 11 3 IMPLICIT INPUTS: 
0000 113; NONE 
0000 114; 
44 12 ; OUTPUT PARAMETERS: 
00000004 0000 119 : result = 4 ; The G COMPLEX*16 result, passed by 
0000 118 3; reference. 
0900 119; 
0000 120 ; IMPLICIT OUTPUTS: 
44 121; NONE 
000 1 § 3 
44 1 3; COMPLETION CODES: 
000 124; NONE 
44 125; 
000 : § : SIDE EFFECTS: 
0000 : 8 ; SS$_ROPRAND If either part of argument is reserved operand. 
3H 
000 132 
007C 0 1 § -ENTRY MTHSCGSQRT, “M<R2, R3, R4, RS, ROD 
; 134 MTHSFLAG_JACKET ; flag as math routine 
6D 00000000'GF 9E MOVAB G*MTHSSJACKET_HND, (FP) 
3; set handler address to jacket 
2 ; handler 
9 135 
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@arg(AP), R2 
#*x8000, R2 


rg 


oa 
> 2 
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:19:1¢ par /yes Macro VO 


arg(AP 
#1, G*MTHSCGABS 
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| 

act zr 

R2-R5 = ABS(r) 

px address of arg on stack 

RO-R1 = CABS((r, 7)) 

RO-R1 = ABS(r) + CABS((r, 7)) 
= (ABS(r) + CABS((r, 1))) / 2 

“Ri = ROOT = SORT (above) 

no, go ahead 

make zero quotient 


e 
= i / ROOT 

=Q= i/ (2 * ROOT) 
r positive 
en return (ROOT, Q) 
se switch ROOT and Q 
i positive 

n return (Q, ROOT) 
else negate ROOT anG Q 
and return (-Q, -ROOT) 


continue to ep ROOT and Q 
and return (Q, ROOT) 


result address 
real part 
imaginary part 


MTHSCGSQRT "9-3 
Symbol table =§ 


ARG = 00000008 
MTHSSJACKET_HND teeeeere 
MTHSCGABS etereree x 
MTHSCGSORT 00000000 RG 1 
MTHSGSORT_RS eereeree xX 
RCSULT = 9000004 

RETRN OOOO0SF R 0 

0000059 R 
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treme ene en amon aae + 

! Psect synopsis ! 

$¢eeecewa sear oceceaaan + 
PSECT name Allocation PSECT No. Attributes 
- ABS. 00000000 0.) 00 ¢ 0.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
_MTHSCODE QOOO0006A ¢ 106.) O01 ¢ 1.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 


E 
Initialization 37 00:00:00.08 00:00:00.58 
Command processing 122 00:00:00.64 BF DR ge 
Pass 1 89 00:00:00.72 Sb STS 7 
Symbol table sort 00:00:00.00 Sb SF -00 
Pass 2 45 00:00:00.50 BOP Rede 
Symbol table output Ses Be 4 00:00:00.23 
Psect synopsis output 00:00:00.0 00:00:00.02 
Cross-reference output 00:00:00.0 00:00:00.00 
Assembler run totals 299 00:00:01.98 00:00:09.81 


The working set Limit was 900 pages. 

2877 bytes (6 pages) of virtual memory were used to buffer the intermediate code. 

There were 10 pages of symbol table space allocated to hold 8 non-local and 2 local symbols. 
9 source Lines were read in Pass 1, producing 11 object records in Pass 2. 

1 page of virtual memory was used to define 1 macro. 
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; Macro Library statistics H 


pe reg werner 
(uacccscncesesencececeese 
Phase Page faults CPU Time Lapsed Time 
Macro Library name Macros defined 
“S255SDUAZB:(SYSLIBISTARLET.MLB;20 0 tS 0 | 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 
MACRO/ENABLE=SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHCGSQRT/OBJ=OBJ$:MTHCGSQRT MSRC$:MTHJACKET/UPDATE=(ENHS: MTHJACKET) +S 
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